Shear viscosity of the quark matter 
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We discuss shear viscosity of the quark matter by using Kubo formula. The shear viscosity is 
calculated in the framework of the quasi-particle RPA for the Nambu-Jona-Lasinio model. We 
obtain a formula that the shear viscosity is expressed by the quadratic form of the quark spectral 
function in the chiral symmetric phase. The magnitude of the shear viscosity is discussed assuming 
the Breit-Wigner type for the spectral function. 
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I. INTRODUCTION 



The existence of the quark gluon plasma (QGP), which 
is predicted by Quantum chromodynamics, has not been 
discovered in Nature. In order to produce such a new 
state of matter, the experimental investigation started at 
the Relativistic Heavy Ion Collider (RHIC). The data at 
RHIC, however, seems to reveal some unexpected prop- 
erties of the high density matter produced in the exper- 
iment [J- [II]. Namely it could be explained by a fluid 
model with small viscosity; it is almost perfect fluid. 
This fact encourages many researchers to calculate the 
transport coefficients of the quark matter [f|-[2(|. How- 
ever their calculations are so complex that the numerical 
results are not settled completely. On the other hand, 
man y st udies by the lattice QCD have been in progress 
[2ll|-|22l] Recently an approach from a very quite differ- 
ent field, black hole physics, has generated a great deal 
of interest. This theory is based on the AdS/CFT cor- 
respondence and predicts the lower bound of the shear 
viscosity of the quark matter (23j-[25|. 



It is the purpose of this paper to calculate the shear vis- 
cosity of the quark matter using the Nambu-Jona-Lasinio 
(NJL) model. It is known that free quark gas has infi- 
nite shear viscosity and the strong interaction between 
quarks lowers the value. We are interested in the physi- 
cal origin of the small viscosity of a strongly interacting 
system. Hence we take up only the quark sector in the 
QGP in this paper. As for the calculation method of the 
viscosity, we use the correlation function method which 
is so-called Kubo formula [26|-[29(. 

In the next section, the Kubo formula for the shear 
viscosity is reviewed. Then it is calculated with the use 
of the quasi-particle RPA in §3. Finally assuming the 
Breit-Wigner type for the spectral function, the shear 
viscosity is evaluated and several discussions are given. 



II. KUBO FORMULA 
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According to the Kubo formula, the shear viscosity 
rj{pS) at temperature T is given by 

V(U) = ^J Q dte 1 ^ J dr(J xy (r,t), J xy (0,0)), (1) 

where J xy is the x, y component of the energy- momentum 
tensor of the quark matter. The correlation function in 
the right-hand side is defined by 



(A,B) = \ d\(e XH Ae- XH B), 
Jo 



(2) 



where A and B are operators of any physical quan- 
tity and H denotes our Hamiltonian. The bracket, 
(A) = Tx(.Ae - ^ )/Tre - ^ , means the thermal average 
at temperature T (/3 = 1 /T) . Using partial integration in 
the right-hand side of Eq.(l), the viscosity can be trans- 
formed into 

77M = -[n R H-n R (o)]. (3) 

Here U r (cj) is a retarded Green's function defined by 
n R {Lo) = -i( dte^ 4 f dr([J xy (r,t), J xy (0,0)}), (4) 



where [ , ] in the integrand denotes the commutation 
relation. Noting that (U r (uj))* = U R (-ui), the (static) 
viscosity is reduced to 



r] = r)(uj = 0) = — -imn fl (w; 

duj 



(5) 



LU = + 



In order to calculate the above H R (ui), it is convenient to 
transform into the imaginary time (Matsubara) formal- 
ism. We introduce the following correlation function, 



n(ja;„) = — I dre 
Jo 
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dr(T T (J xy (r,T)J xy (0,0))), 
(6) 

where the Matsubara frequency is represented by Lu n = 
2irnT and T T means the (imaginary) time ordering op- 
erator. As is well known, the retarded Green's func- 
tion H r (uj) is obtained by the analytic continuation: 





FIG. 1: The diagram for the interaction with T = 1 or i 7 5T 

We take the Nambu-Jona-Lasinio (NJL) model for the 
quark matter in this paper [30( • The Lagrangian density 
is given by 

C = V;(i7 • d - m)i> + g[(V>V) 2 + (V^75tV>) 2 ]- (7) 

Here is the field operator for quarks and the quark 
mass is assumed to be zero (m = 0) . Then the canonical 
energy-momentum tensor is read as 

J*, = 5[W-W# (8) 

If this expression is substituted into Eq.(6), the corre- 
lation function LT is calculated according to the Wick 
theorem (Feynman diagram method). Our interaction 
is represented by the Feynman diagram Fig.l where T is 
1 (unit matrix) or ijsT. Note that the four Fermi oper- 
ators are composed of two pairs which are connected by 
a broken line. 



III. QUASI-PARTICLE RPA 

Here we take the Fermi liquid theory and the ran- 
dom phase approximation (quasi-particle RPA) . Namely 
the correlation function is approximated by the ring di- 
agrams as shown in Fig. 2. All the propagators must be 
the dressed ones because the free propagator gives rise to 
the infinity in the lowest term as shown in the last sec- 
tion. Therefore the ring approximation is broken down. 
In fact the quarks have strong correlation even in the chi- 
ral symmetric phase [3l| . A similar calculation was done 
in a relativistic scalar field theory in the Ref.(30) where 
the shear viscosity was calculated by taking into account 
chain diagrams. 

The present approximation is also derived by the 1 /N 
expansion used in QCD (N = N c ) Let a dia- 

gram has n vertices and ni quark loops. Let us consider 
the ,/V-dependence of this diagram. Each quark loop gives 
rise to N from the trace of color variable. The coupling 
constant of the NJL model is order of 1/N because the 
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FIG. 2: The ring diagrams for the correlation function II. 

quark-gluon coupling constant is order of 1/yN in the 
1/N expansion. Consequently this diagram is order of 
TV™' - ™. Therefore the leading diagrams are those with 
the maximum number of the loop for the fixed n. It is 
proved that ni < n + 1. 

proof. Let us consider any diagram with n vertices. If the 
number of the quark propagators is rij, we have the rela- 
tion, 2rii — 4 + An. Each propagator belongs to only one 
loop necessarily according to the Wick theorem. Since 
each loop contains two or more propagators, the maxi- 
mum number of the loop is realized that all the loops are 
composed of two propagators: n; = rii/1 = n + 1. The 
loop composed of one propagator (tadpole) is included 
in the dressed propagator. Thus it is evident that such 
diagrams correspond to the ring ones of Fig. 2 and their 
mag nitude is 0(N V ). 

For example, let us consider the first order diagrams 
with respect to the interaction shown in Fig. 3. One has 
two loops (Hartree term) and the other has only one loop 
(Fock term) composed of four propagators. The former 
(a) is 0(N 1 ) because it has one vertex and two quark 
loops. The latter (b) with one vertex and one loop is 
0(1). Consequently the latter term should be neglected. 
In the same way the third diagram in Fig. 2 is order of N 
and so on. 

Now the correlation function is written as 

n (^») = \Y.j ^3^ 2 Tr[ 7 2 G(p,^+c„) 7 2 

G( P , WI )] + -^y —— Mx 

Tr[ 7 2 G(p,^+ W „)rG(p,^)] 

Tr[rG(q, LJ m - cu„) 7 2 G(q, u m )] • • • , (9) 

where G(p,iuj n ) is the dressed quark propagator. It 
should be noted that the dressed propagator is propor- 
tional to the linear combination of 7 M in the chiral limit 
;36]. This leads to Tr[ 7 2 GrG] = in the second term be- 
cause the trace has an odd product of 7 matrices (r = 1 
or i 7 5r). Thus the second term vanishes. Similarly the 
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z= e 



FIG. 3: The first order diagrams: (a) Hartree term and (b) 
Fock term. 



FIG. 4: The contour for the calculating the integral of z in 
Eq.(ll). 



third term contains the same loop with odd product of 
7 matrices so that it vanishes. The same calculation is 
realized in all the higher order terms in Fig. 2. Conse- 
quently the correlation function II is reduced to only the 
first term in our present case. 

In order to calculate the correlation function, let us 
follow the procedure taken for the calculation of the elec- 
trical conductivity in Ref[35]. We consider the spectral 
representation for the dressed propagator, which is writ- 
ten as 



de p a/3 (p,e) 



(10) 



We substitute this expression into the correlation func- 
tion (Eq.(8)) and the summation over Matsubara fre- 
quency is replaced by the contour integral: 

S = T^Tr[ 7 2 G(p,^ / +^„) 7 2 G(p, i ^)] 
i 

f dz 

= - n(z)Tr[G(p,z) 7 2 G(p,z + i ^)7 2 ](H) 

where n(z) = (1 + e l3z )~ 1 is the Fermi distribution func- 
tion. The contour G is divided into three pieces as shown 
in Fig. 4, because the integrand has branch cuts on the 
two lines z — e and z — e — iu> n where e is real (Note 
that poles of n(z) do not lie on the two branch cuts). 

Since the integral along the large circle vanishes, it is 
rewritten as 



S = - 



d£ 

2ni 



T n(£)Tr[G(e + i6)j 2 G(e + iuj n )-f 2 



-G(e - i6)j 2 G(e + ioj n h 2 + G(e - iw„)7 2 
G(e + iS)j 2 - G(e - iuj n )j 2 G(e - ^) 7 2 ], 

where S is an infinitesimal positive number introduced 
in order to avoid the branch cuts. Here the argument 
p is abbreviated for convenience. Noting the relation, 
G(e + iS) — G(e — id) = —ip(e), the above equation is 



expressed as follows: 

°° d£ 



S 



2n 



n( £ )Tr[(G(e + iw n ) + G(e - iuj n ))Y p{e)Y 



(12) 

Substituting the spectral representation (10) again, wc 
get 



S = 



^£ r — 

^ J-oo 2^ 

Tr[p(e'h 2 p(eh 2 } 



1 



1 



e + iu>„ — e' 



e — ioj„ — e 

(13) 



By the analytic continuation iu> n — > u) + iS, the imagi- 
nary part of the above equation becomes 



ImS = 



~{n(e + u:)-n{e))TMe + ^)l 2 p(e)l\ 

Z7T Z 



(14) 

because Tr[p(e')7 2 /9(e)7 2 ] is real. As a result we obtain 
an expression for the shear viscosity which is related to 
the quark spectral function, 



1 



V = 



d£ 

2vr 



dp 

(27T) 



,dn r 



3 P*~— Tr [p(P> e )7 p{p,e)l 



de 



(15) 

This equation means that the calculation of the shear 
viscosity is reduced to that of the spectral function. The 
factor of dn/de means that the main contribution comes 
from the neighborhood of the Fermi surface. Note that 
the information of the Fermi surface (chemical potential) 
is contained in the quark spectral function. 



IV. NUMERICAL RESULTS 

Now we are in a position to discuss the quark spectral 
function when the chiral symmetry is restored. Accord- 
ing to Weldon [3(| , the general forms of the self-energy 
for the massless quark is expressed by 



S(p) = S +7 °A + + I]_7 A_, 



(16) 
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FIG. 5: The shear viscosity as a function of the imaginary 
part of the self-energy at T = 200MeV and /i = lOMeV 



using the projection operator A± = (1 ± pr) and po = 
e + [i (fi: the chemical potential). Then the retarded 
(advanced) Green's function are given by 



G R (p) = 
G A (p) = 



7° A H 



7° A. 



P°+ - P - S. 
■+ 



P°+ 



7° A H 



■p-E H 
7 °A_ 



P°+-P-xt pI+p-^x 

The spectral function is given by 

pip) = i[G R G A ] = P+ 7°A+ + p_7°A_, 



(17) 



(18) 



FIG. 6: The ratio of the shear viscosity to the entropy as 
a function of the imaginary part of the self-energy at T — 
200MeV and fj, = lOMeV 



soft mode 



where the p± is defined by 

P±{P) " (19) 

Thus we get a spectral function of the Breit-Wigner type. 

As for the calculation of the self-ener gy S , the coupling 
to the soft mode [13] or the plasmino |38T|-[4l| might be 
probable. Before that it is instructive to calculate the 
shear viscosity by using a simple form, 

E±(p) = M-*T, (20) 

where M and T are positive constants and regarded as 
free parameters. Since the quasi-particles near the Fermi 
surface contribute main part of the integral in Eq.(15), 
this replacement is not so bad. The real part M is taken 
to be zero since it essentially shifts the chemical potential. 
The imaginary part T means the inverse of the mean free 
path (or the inverse of the life time) of the quasi-particle. 
The numerical calculation was done with the temperature 
T = 200MeV and the chemical potential n = lOMeV for 
convenience. The shear viscosity is shown in Fig. 5 as a 
function of T and the ratio of the shear viscosity to the 
entropy density is also shown in Fig. 6. 

These figures evidently show that the viscosity is a 
rapidly increasing function of the mean free path, 
The viscosity seems to diverge at T = where the quark 



FIG. 7: The self-energy diagram of the quark which is caused 
by the coupling to the soft modes. 



matter becomes free gas (the mean free path is infinite). 
This is understood by noting that (T/(x 2 + T 2 )) 2 — > 
(tt5(x)) 2 as r — > 0. As the T becomes larger, the vis- 
cosity decreases rapidly. This is also apparent because 
(r/(x 2 + T 2 )) 2 ->■ as r -> oo. This V dependence of 
the viscosity is consistent with the classical expression of 
the viscosity of gas: 77 = pvl/3 (p=density, u=velocity 
and 1= mean free path of the particle). The ratio rj/s 
of water in the normal state is about 30. Therefore the 
quark matter is a perfect fluid under strongly interact- 
ing state (r > lMeV). The lower bound predicted by 
the superstring theory is drawn by the dashed line in the 
Fig.6. 

Next let us comment the evaluation of the T by consid- 
ering the coupling between the quark and the soft mode 
as shown in Fig. 7 [3fJ,[? ]. The dashed line denotes 
for the propagator of the soft mode. As an example, we 
show here T(p = 0, w) shown in Fig.8 at T = 200MeV and 
fx = lOMeV. (The critical temperature is T c = 165MeV) 
This figure suggests that the quark matter has strong 
correlation even in the chiral symmetric phase: the T is 
of order 0.1 GeV. According to Fig. 5, the shear viscosity 
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of the quark matter with the use of the Kubo formula. 
Assuming the quasi-particle RPA, we have obtained a 
simple formula for the shear viscosity, which is expressed 
by the quadratic form of the quark spectral function in 
the chiral symmetric phase. The magnitude of the shear 
viscosity is expected to be small in the strongly interact- 
ing system. 
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FIG. 8: The calculated V(p = 0,w) at T = 170MeV and 
H = lOMeV 

would be small enough. Of course the detailed calcula- 
tion of the shear viscosity itself is left in the future. 
In conclusion, we have calculated the shear viscosity 
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